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Traditionally, studies aimed at inferring the distribution of birth periods of neu- 
tron stars are based on radio surveys. Here we propose an independent method to 
constrain the pulsar spin periods at birth based on their X-ray luminosities. In par- 
ticular, the observed luminosity distribution of supernovae poses a constraint on the 
initial rotational energy of the embedded pulsars, via the Lx — E rot correlation found 
for radio pulsars, and under the assumption that this relation continues to hold beyond 
the observed range. We have extracted X-ray luminosities (or limits) for a large sample 
of historical SNe observed with Chandra, XMM and Swift, that have been firmly clas- 
sified as core-collapse supernovae. We have then compared these observational limits 
with the results of Monte Carlo simulations of the pulsar X-ray luminosity distribu- 
tion, for a range of values of the birth parameters. We find that a pulsar population 
dominated by millisecond periods at birth is ruled out by the data. 



1 INTRODUCTION 



Modeling the observed properties of the Galactic population 
of radio pulsars, with the purpose of inferring their intrinsic 
properties, has been the subject of extensive investigation 
for several decades (e.g. Gunn & Ostriker 1970; Phinney & 
Blandford 1981; Lyne et al. 1985; Stollman 1987; Emmer- 
ing & Chevalier 1989; Narayan & Ostriker 1990; Lorimer et 
al. 1993; Hartman et al. 1997; Cordes & Chernoff 1998; Ar- 
zoumanian, Cordes & Chernoff 2002; Vranesevic et al 2004; 
Faucher-Giguere & Kaspi 2006; Ferrario & Wickramasinghe 
2006). Since the fraction of pulsars that can be detected 
close to their birth constitutes a negligible fraction of the 
total sample, these studies generally use the present day ob- 
served properties of pulsars (namely their period P and pe- 
riod derivative P), together with some assumptions about 
their time evolution, to reconstruct the birth distribution 
of periods and magnetic fields for the pulsar population. 
These analyses also need to make assumptions about pul- 
sar properties and their evolution (such as, for example, 
the exact shape of the radio beam and its dependence on 
the period), as well as overcome a number of selection ef- 
fects. Results from various investigations have often been 
conflicting, with some studies favoring initial periods in the 
millisecond range (e.g. Arzoumanian et al. 2002), and oth- 
ers instead finding more likely periods in the range of sev- 
eral tens to several hundreds of milliseconds (e.g. Faucher- 



Giguere & Kaspi 2006). The efforts put over the years into 
this area of research stem from the fact that the birth prop- 
erties of neutron stars (NSs) are intimately related to the 
physical processes occurring during the supernova (SN) ex- 
plosion and in the proto-neutron star. As such, they bear 
crucial information on the physics of core-collapse SNe, in 
which most are thought to be formed. 

Besides the inferences on the birth parameters from the 
radio population discussed above, we show here that con- 
straints can be derived also from the X-rays. Young, fast 
rotating neutron stars are indeed expected to be very bright 
in the X-rays. In fact, observationally there appears to be 
a correlation between the rotational energy loss of the star, 
-Erot, and its X-ray luminosity, L x . This correlation was no- 
ticed by Verbunt et al. (1996), Becker & Trumper (1997), 
Seward & Wang (1988), Saito (1998) for a small sample of 
objects, and later studied by Possenti et al. (2002; P02 in the 
following) for the largest sample of pulsars known to date. 

Combining the birth parameters derived from the ra- 
dio (which determine the birth distribution of E Iot for the 
pulsars), with the empirical L x — E YOt correlation, the dis- 
tribution of X-ray luminosity can be predicted for a sam- 
ple of pulsars with a certain age distribution. The above 
calculation was performed by Perna & Stella (2004). They 
found that the birth parameters derived by Arzoumanian et 
al. (2002), together with the L x — E TOt correlation derived 
by P02, yield a sizable fraction of sources with luminosities 
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> 10 39 erg/s, which could hence constitute potential contrib- 
utors to the observed population of ultra luminous X-ray 
sources (ULXs) observed in nearby galaxies (e.g. Fabbiano 
& White 2003; Ptak & Colbert 2004). Obviously, these pre- 
dictions were heavily dependent on the assumed initial birth 
parameters (the periods especially) of the pulsar population. 

In this paper, we propose a new, independent method 
to constrain the pulsar spin periods at birth from X-ray ob- 
servations, and hence also assess the contribution of young, 
fast rotating NSs to the population of bright X-ray sources. 
Since neutron stars are born in supernova explosions, and 
very young pulsars are still embedded in their supernovae, 
the X-ray luminosity of the SNe provides an upper limit 
to the luminosity of the embedded pulsars. We have ana- 
lyzed an extensive sample of historical SNe whose position 
has been observed by Chandra, XMM or Swift, and studied 
their X-ray counterparts. We measured their X-ray lumi- 
nosities, or derived a limit on them in the cases of no detec- 
tion. A comparison between these limits and the theoretical 
predictions for the distribution of pulsar X-ray luminosities 
shows that, if the assumed initial spins are in the millisecond 
range, the predicted distribution of pulsar X-ray luminosi- 
ties via the L x — _Z5 ro t correlation is highly inconsistent with 
the SN data. Our analysis hence suggests that a substantial 
fraction of pulsars cannot be born with millisecond periods. 

The paper is organized as follows: in §2, we describe the 
method by which the SN X-ray flux measurements and limits 
are extracted, while in §3 we describe the theoretical model 
for the distribution of the X-ray luminosity of young pulsars. 
A comparison between the theoretical predictions and the 
data is performed in §4, while the results are summarized 
and discussed in §5. 



2 X-RAY ANALYSIS OF HISTORICAL 

SUPERNOVAE OBSERVED BY CHANDRA, 
XMM AND SWIFT 

We compared and combined the CfA List of SupernovaqH 
the Padova-Asiago Catalogued, the Sternberg Catalogue^ 
(Tsvetkov et al. 2004), and Michael Richmond's Super- 
nova PagcQ, to create a list of unambiguosly identified core- 
collapse SNe (updated to 2007 April). We cross-correlated 
the SN positions with the catalogues of Chandra/ ACIS, 
XMM-Newton/EPIC and Swift/XRT observation^, to de- 
termine which SN fields have been observed by recent X-ray 
missions (ASCA was excluded because of its low spatial res- 
olution, and ROSAT because of its lack of 2-10 keV sensi- 
tivity) . For the Chandra ACIS-S data, we limited our search 
to the S3 chip. We obtained a list of ~ 200 core-collapse SNe 
whose positions happened to be in a field observed at least 
once after the event. From the list, we then selected for this 
paper all the core collapse SNe with unambiguos subtype 
classification (Type Ib/c, Type Iln, IIL, and IIP and lib). 

1 http://cfa-www.harvard.edu/iau/lists/Supernovae.html, com- 
piled by The Central Bureau for Astronomical Telegrams at the 
Harvard-Smithsonian Center for Astrophysics. 

2 http:/ /web. pd.astro.it/supern/snean.txt 

3 VizieR On-line Data Catalog: 11/256 

4 http:/ /stupendous. rit.edu/richmond/sne/sn.list 

5 Search form at http://heasarc.nasa.gov 



That is about half of the total sample. We leave the analysis 
of the other ~ 100 SNe (classified generically as Type II) to 
a follow-up paper. 

We retrieved the relevant X-ray datasets from the pub- 
lic archives of those three missions. The optical position of 
each SN in our sample is well known, to better than 1": 
this makes it easier to determine whether a SN is detected 
in the X-ray band (in particular for Chandra), even with a 
very low number of counts, at a level that would not be con- 
sidered significant for source detection in a blind search. For 
the Chandra observations, we applied standard data analysis 
routines within the Chandra Interactive Analysis of Obser- 
vations (CIAO) software packag^fl version 3.4. Starting from 
the level- 2 event files, we defined a source region (radius 27 5, 
comprising ~ 95% of the source counts at 2 keV, on axis, and 
proportionally larger extraction radii for off-axis sources) 
and suitable background regions not contaminated by other 
sources and at similar distances from the host galaxy's nu- 
cleus. For each SN, we extracted source and background 
counts in the 0.3-8 keV band with dmextract. In most cases, 
we are dealing with a very small number of counts (e.g., 2 or 
3, inside the source extraction region) and there is no excess 
of counts at the position of the SN with respect to the lo- 
cal background. In these cases, we calculated the 90% upper 
limit to the number of net counts with the Bayesian method 
of Kraft et al. (1991). We then converted this net count-rate 
upper limit to a flux upper limit with WebPimmfl, assum- 
ing a power-law spectral model with photon index F = 2 
and line-of-sight Galactic column density. The choice of a 
power-law spectral model is motivated by our search for X- 
ray emission from an underlying pulsars rather than from 
the SN shock wave. In a few cases, there is a small excess of 
counts at the SN position: we then also built response and 
auxiliary response functions (applying psextract in CIAO), 
and used them to estimate a flux, assuming the same spec- 
tral model. When possible, for sources with ~ 20-100 net 
counts, we determined the count rates separately in the soft 
(0.3-1 keV), medium (1-2 keV) and hard (2-8 keV) bands, 
and used the hard-band rates (essentially uncontaminated 
by soft thermal-plasma emission, and unaffected by the un- 
certainty in the column density and by the degradation of 
the ACIS-S sensitivity) alone to obtain a more stringent 
value or upper limit to the non-thermal power-law emission. 
Very few sources have enough counts for a two-component 
spectral fit (mekal thermal plasma plus power-law): in those 
cases, we used the 2-10 keV flux from the power-law com- 
ponent alone in the best-fitting spectral model. For those 
spectral fits, we used the XSPEC version 12 software pack- 
age (Arnaud 1996). 

When we had to rely on XMM-Newton/EPIC data, we 
followed essentially the same scheme: we estimated source 
and background count rates (this time, using a source ex- 
traction circle with a 20" radius) in the full EPIC pn and 
MOS bands (0.3-12 keV) and, when possible, directly in the 
2-10 keV band. The count rate to flux conversion was ob- 
tained with WebPimms (with a F — 2 power-law model ab- 
sorbed by line-of-sight column density) or through full spec- 



6 http://cxc.harvard.edu/ciao 

7 http://heasarc.nasa.gov/Tools/w3pimms.html and 
http: / /cxc. harvard.edu/toolkit/pimms.jsp 
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SN Host galaxy Type Age (yr) £2-10 fcev( er g/ s ) Instrument Observation date 



1923A 


N5236 


IIP 


77.3 


< 6.0 x 10 35 


1926A 


N4303 


IIP 


75.3 


< 1.4 x 10 37 


1937A 


N4157 


IIP 


67.3 


< 1.3 x 10 37 


1937F 


N3184 


IIP 


62.1 


< 2.7 x 10 36 


1940 A 


N5907 


ILL 


63.0 


< 1.0 x 10 37 


1940B 


N4725 


IIP 


62.6 


< 8.6 x 10 36 


1941A 


N4559 


IIL 


60.2 


< 5.5 x 10 36 


1948B 


N6946 


IIP 


55.1 


< 4.7 x 10 35 


1954A 


N4214 


lb 


48.9 


< 1.6 x 10 35 


1959D 


N7331 


IIL 


41.6 


< 2.2 x 10 37 


1961V 


N1058 


I In 


38.3 


< 6.1 x 10 37 


1962L 


N1073 


Ic 


41.2 


< 4.7 x 10 37 


1962M 


N1313 


IIP 


40.3 


< 3.7 x 10 3S 


1965H 


N4666 


IIP 


37.7 


< 1.5 x 10 38 


1965L 


N3631 


IIP 


37.8 


< 5.7 x 10 36 


1968L 


N5236 


IIP 


32.0 


< 1.5 x 10 36 


1969B 


N3556 


IIP 


32.6 


< 3.8 X W ou 


1969L 


N1058 


IIP 


30.3 


< 4.8 X 10 


1970G 


N5457 


IIL 


33.9 


< 4.9 X 1U 


1972Q 


TVT A r» r A 

N4254 


IIP 


30.5 


< 3.0 x 10 


1972R 


"MOO A "1 

N2841 


lb 


31.9 


< 7.3 x 10 


1973R 


N3627 


IIP 


25.9 
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< 7.7 X 10 


1976B 


AT A A C\i\\ 
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< 8.9 X 10 


1979C 
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1980K 
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< 6.5 x 10 36 
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1 n n 
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IC 
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^ 1 1 v 1 n36 


lyooiN 


IN OZoD 


TK 
ID 


1 ft fi 
lO.o 
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1985L 




IIL 


14.9 


< 8.1 x 10 37 


1986E 


N4302 


IIL 


19.6 


1 4 + S 5 x 10 38 


19861 


N4254 


IIP 


17.1 


<• 5(1 y 1 A 38 


1986J 


N891 


I In 


21.2 


8 5+ ' 5 v 10.38 

O.O q r A 1U 


1986L 


N1559 


IIL 


18.9 


< 1.4 X 10 38 


1987B 


N5850 


Iln 


14.1 


< 1.5 x 10 38 


1988A 


N4579 


IIP 


12.3 


< 2.4 x 10 37 


1988Z 


MCG+03-28-22 


Iln 


15.5 


2 9 +0 - 5 X 10 39 




in / 4 1 y 


Ic 
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1 1 + -k v m 39 

J-'-L — O 5 x 1U 


1991N 
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1996ac 
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< 6.1 x 10 37 


1 996hu 


N3631 


Iln 


6.6 


< 2.1 x 10 37 


1996cr 


ES097-G13 


Iln 


4.2 


1 9 +0 ' 4 x 10 39 

— 0.4 


1997X 


N4691 


Ic 


6.1 


< 2.2 x 10 37 


1997bs 


N3627 


Iln 


2.5 


< 2.9 x 10 38 


1998S 


N3877 


Iln 


3.6 


3.8+!^ x 10 39 

— 0.5 


1998T 


N3690 


lb 


5.2 


<- 2 n x in 38 


1 yyo uw 
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0.0 
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TTvi 
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< O.D X 1U 


1 999om 


N1637 


IIP 


1.0 


< 1.4 x 10 37 


1999gi 


N3184 


IIP 


0.10 


2.6t°;« x 10 37 


2000P 


N4965 


Iln 


7.2 


< 1.2 x 10 39 


2000bg 


N6240 


Iln 


1.3 


< 1.4 x 10 39 


2001ci 


N3079 


Ic 


2.5 


< 5.0 x 10 37 


2001du 


N1365 


IIP 


1.3 


< 3.8 x 10 37 


2001em 


UGC11794 


Ib/Ic 


4.7 


5.8±i"| x 10 40 


2001gd 


N5033 


lib 


1.1 


1.0+° :3 x 10 39 


2001ig 


N7424 


lib 


0.50 


3.5+^ x 10 37 



ACIS 2000-04-29, 2001-09-04 

ACIS 2001-08-07 

EPIC 2004-05-16 

ACIS 2000-01-08, 2000-02-03 

EPIC 2003-02-20, 2003-02-28 

ACIS 2002-12-02 

ACIS 2001-01-14, 2001-06-04, 2002-03-14 

ACIS 2001-09-07, 2002-11-25, 2004-10-22, 2004-11-06, 2004-12-03 

ACIS 2003-03-09 

ACIS 2001-01-27 

ACIS 2000-03-20 

ACIS 2004-02-09 

ACIS 2002-10-13, 2002-11-09, 2003-10-02, 2004-02-22 

ACIS 2003-02-14 

ACIS 2003-07-05 

ACIS 2000-04-29, 2001-09-04 

ACIS 2001-09-08 

ACIS 2000-03-20 

ACIS 2004-07-05, 2004-07-11 

EPIC 2003-06-29 

EPIC 2004-11-09 

ACIS 1999-11-03 

EPIC 2002-07-01 

ACIS 2006-02-18 

ACIS 2004-10-22, 2004-11-06, 2004-12-03 

ACIS 2004-07-29, 2004-11-20 

EPIC 2001-11-24, 2002-05-10 

ACIS 2001-02-06 

ACIS 2000-04-29 

ACIS 2002-12-24 

ACIS 2000-04-28 

ACIS 2005-12-05 

EPIC 2003-06-29 

ACIS 2003-12-10 

EPIC 2005-08-10, 2005-10-12 

EPIC 2001-01-25, 2001-08-26 

ACIS 2000-05-02 

ACIS 2004-06-29 

EPIC 2001-06-19 

ACIS 2003-01-25 

ACIS 2001-04-22 

ACIS 2000-06-20, 2001-06-23, 2003-08-07 

ACIS 2002-05-17 

ACIS 2004-03-27 

ACIS 2002-04-05 

ACIS 2003-07-05 

ACIS 2000-03-14 

ACIS 2003-03-08 

ACIS 1999-11-03 

ACIS 2001-10-17 

ACIS 2003-04-30 

ACIS 2001-10-27 

ACIS 2004-01-25 

EPIC 2005-08-31 

EPIC 2005-04-30, 2005-06-05 

ACIS 2000-10-30 

ACIS 2000-01-08, 2000-02-03 

XRT 2007-05-16 

ACIS 2001-07-29 

EPIC 2003-10-14 

ACIS 2002-12-24 

EPIC 2006-06-14 

EPIC 2002-12-18 

ACIS 2002-06-11 
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SN Host galaxy Type Age (yr) £2-10 keV ( cr g/ s ) Instrument Observation date 



2002ap 


N628 


Ic 


0.92 


< 3.1 x 10 36 


EPIC 


2003-01-07 


2002fj 


N2642 


Iln 


4.7 


< 1.3 x 10 39 


XRT 


2007-05-11 


2002hf 


MCG-05-3-20 


Ic 


3.1 


< 7.6 x 10 3S 


EPIC 


2005-12-19 


2003L 


N3506 


Ic 


0.08 


T.7 + _\\ x 10 39 


ACIS 


2003-02-10 


2003ao 


N2993 


IIP 


0.016 


5.6+^ x 10 38 


ACIS 


2003-02-16 


2003bg 


MCG-05-10-15 


Ic/IIb 


0.33 


x 10 38 


ACIS 


2003-06-22 


2003dh 


Anon. 


Ic 


0.71 


< 5.0 x 10 4() 


EPIC 


2003-12-12 


2003jd 


MCG-01-59-21 


Ic 


0.041 


< 3.0 x 10 3S 


ACIS 


2003-11-10 


20031w 


Anon. 


Ic 


0.34 


7.0l 3 x 10 40 


ACIS 


2004-04-18 


2004C 


N3683 


Ic 


3.1 


1.0+ 9 ' 3 x 10 38 


ACIS 


2007-01-31 


2004dj 


N2403 


IIP 


0.34 


l.l+° :3 x 10 37 


ACIS 


2004-12-22 


2004dk 


N6118 


lb 


0.030 


1.5+ ; 4 x 10 39 


EPIC 


2004-08-12 




N6946 


IIP 


0.18 


1 +0 - 2 x 10 38 


ACIS 


2004-10-22, 2004-11-06, 2004-12-03 




N5420 


Th fir 


U. JU 


s 1 n y 1 f|40 

<^ l.U X _LU 


XRT 


onn^ n7 1 7 

iUUO-U I -1 I 






lib 


0.041 


< 1.1 x 10 39 


ACIS 


2005-02-1 4 


2005at 


N6744 


Ic 


1.7 


< 3 X 1 3S 


XRT 


2006-10-31 


2005bf 


MCG+00-27-5 


lb 


0.58 


< 6.0 x 10 39 


XRT 


2005-11-07 


2005bx 


MCG+12-13-19 


Iln 


0.25 


< 1.0 x 10 39 


ACIS 


2005-07-30 


2005da 


UGC11301 


Ic 


0.098 


< 5.0 x 10 39 


XRT 


2005-08-23 


2005db 


N214 


Iln 


0.036 


< 2.0 x 10 39 


EPIC 


2005-08-01 


2005ck 


UGC2526 


Ic 


0.041 


< 4.0 x 10 39 


XRT 


2005-10-07 


2005gl 


N266 


Iln 


1.6 


< 3.4 x 10 39 


XRT 


2007-06-01 


2005kd 


Anon. 


Iln 


1.2 


2.6l°; 4 x 10 41 


ACIS 


2007-01-24 


2006T 


N3054 


lib 


0.0082 


< 6.0 x 10 39 


XRT 


2006-02-02 


2006aj 


Anon. 


Ic 


0.43 


< 7.0 x 10 39 


XRT 


2006-07-25 


2006bp 


N3953 


IIP 


0.058 


i.o+°i x 10 38 

w — 0.2 


EPIC 


2006-04-30 


2006bv 


UGC7848 


Iln 


0.0082 


< 1.2 x 10 39 


XRT 


2006-05-01 


2006dn 


UGC12188 


lb 


0.033 


< 2.5 x 10 4() 


XRT 


2006-07-17 


2006gy 


N1260 


Iln 


0.16 


< 2.0 x 10 3S 


ACIS 


2006-11-15 


2006jc 


UGC4904 


lb 


0.068 


2.1±° 6 6 x 10 38 


ACIS 


2006-11-04 


20061c 


N7364 


Ib/Ic 


0.016 


< 2.0 x 10 40 


XRT 


2006-10-27 


20061t 


Anon. 


lb 


0.068 


< 4.0 x 10 39 


XRT 


2006-11-05 


2007C 


N4981 


lb 


0.022 


< 3.0 x 10 4() 


XRT 


2007-01-15 


2007D 


UGC2653 


Ic 


0.025 


< 3.0 x 10 4() 


XRT 


2007-01-18 


20071 


Anon. 


Ic 


0.016 


< 9.0 x 10 39 


XRT 


2007-01-20 


2007bb 


UGC3627 


Iln 


0.022 


< 4.4 x 10 39 


XRT 


2007-04-10 



Table 1. X-ray measurements and upper limits for our sample of historical supernovae. When more than one observation was used for 
a given source, the age is an average of three epochs weighted by their exposure lengths. The instrument ACIS is on-board Chandra, 
EPIC on XMM-Newton and XRT on Swift. 



tral analysis for sources with enough counts. We used stan- 
dard data analysis tasks within the Science Analysis Sys- 
tem (SAS) version 7.0.0 (for example, xmmselect for source 
extraction). All three EPIC detectors were properly com- 
bined, both when we estimated count rates, and when we 
did spectral fitting, to increase the signal-to-noise ratio. In 
fact, in almost all cases in which a source position had been 
observed by both Chandra and XMM-Newton, Chandra pro- 
vided a stronger constraint to the flux, because of its much 
narrower point-spread function and lower background noise. 
The Swift data were analyzed using the Swift Software ver- 
sion 2.3 tools and latest calibration products. Source counts 
were extracted from a circular region with an aperture of 20" 
radius centered at the optical positions of the SNe. In some 
cases, Swift observations referred to a Gamma ray burst 
(GRB) associated to a core-collapse SN: we did not obvi- 
ously consider the GRB flux for our population analysis. 
Instead, for those cases, we considered the most recent Swift 
observation after the GRB had faded, and used that to de- 
termine an upper limit to a possible pulsar emission. We 



only considered Swift observations deep enough to detect or 
constrain the residual luminosity to < 10 40 erg s _1 . 

In some cases, two or more Chandra or XMM-Newton 
observations of the same SN target were found in the archive. 
If they were separated by a short interval in time (much 
shorter than the time elapsed from the SN explosion), we 
merged them together, to increase the detection threshold. 
The reason we can do this is that we do not expect the un- 
derlying pulsar luminosity to change significantly between 
those observations. However, when the time elapsed between 
observations was comparable to the age of the SN, we at- 
tributed greater weight to the later observations, for our flux 
estimates. The reason is that the thermal X-ray emission 
from the shocked gas tends to decline more rapidly (over a 
few months or years) than the non-thermal pulsar emission 
(timescale > tens of years). More details about the data 
analysis and the luminosity and color/spectral properties of 
individual SNe in our sample will be presented elsewhere 
(Pooley et al. 2008, in preparation). Here, we are mainly 
interested in a population study to constrain the possible 
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presence and luminosity of high-energy pulsars detectable 
in the 2-10 keV band. 

While this is, to the best of our knowledge, the first X- 
ray search for pulsar wind nebulae (PWNe) in extragalactic 
SNe, and the first work that uses these data to set statisti- 
cal constraints on the properties of the embedded pulsars, it 
should be noted that the possibility of observing pulsars in 
young SNe (a few years old) was originally discussed, from 
a theoretical point of view, by Chevalier & Fransson (1992). 
Furthermore, searches for PWNe in extragalactic SNe have 
been performed in the radio (Reynolds & Fix 1987; Bar- 
tel & Bietenholz 2005). Observationally, however, clear evi- 
dence for pulsar activity in SNe has been lacking. The radio 
emission detected in some SNe, although initially ascribed 
to pulsar activity (Bandiera, Pacini & Salvati 1984), was 
later shown to be well described as the result of circum- 
stellar interaction (Lundqvist & Fransson 1988). There is 
however a notable exception, that is SN 1986J, for which 
the observed temporal decline of the H a luminosity (Rupen 
et al. 1987) has been considered suggestive of a pulsar en- 
ergy input (Chevalier 1987). As noted by Chevalier (1989), 
a possible reason for the apparent low-energy input in some 
cases could be the fact that the embedded neutron stars 
were born with a relatively long period. The present work 
allows us to make a quantitative assessment on the typical 
minimum periods allowed for the bulk of the NS population. 
The list of SNe, their measured fluxes and their ages (at the 
time of observation) are reported in Table 1. 



3 THEORETICAL EXPECTATIONS FOR THE 
X-RAY LUMINOSITY OF YOUNG PULSARS 

Most isolated neutron stars are X-ray emitters throughout 
all their life: at early times, their X-ray luminosity is pow- 
ered by rotation (e.g. Michel 1991; Becker & Trumper 1997); 
after an age of ~ 10 3 — 10 4 yr, when the star has slowed 
down sufficiently, the main X-ray source becomes the inter- 
nal heat of the stai0, and finally, when this is exhausted, 
the only possible source of X-ray luminosity would be ac- 
cretion by the interstellar medium, although to a very low 
luminosity level, especially for the fastest stars (e.g. Blaes 
& Madau 1993; Popov et al. 2000; Perna et al. 2003). An- 
other possible source of X-ray luminosity that has often 
been discussed in the context of NSs is accretion from a 
fallback disk (Colgate 1971; Michel & Dressier 1981; Cheva- 
lier 1989; Yusifov et al. 1995; Chatterjee et al. 2000; Al- 
par 2001; Perna et al. 2000). Under these circumstances, 
accretion would turn off magnetospheric emission, and X- 
ray radiation would be produced as the result of accre- 
tion onto the surface of the star. For a disk to be able 
to interfere with the magnetosphere and accrete, the mag- 
netospheric radius R m ~ 6.6 x 10 7 B^' rn~ 2 ^ 7 cm (with 
m -2//7 being the accretion rate in Eddington units, and 
S12 = B/(10 12 G)) must be smaller than the corotation ra- 
dius R COI ~ 1.5 x 1O 8 P 2/3 (M/M ) 1/3 cm. If, on the other 
hand, the magnetospheric radius resides outside of the coro- 
tation radius, the propeller effect (Illarionov & Sunyaev 



8 In the case of magnetars, this internal heat is provided by mag- 
netic field decay, which dominates over all other energy losses. 



1975) takes over and inhibits the penetration of material 
inside the magnetosphere, and accretion is (at least largely) 
suppressed. For a typical pulsar magnetic field B\2 ~ 5, the 
magnetospheric radius becomes comparable to the corota- 
tion radius for a period P ~ 1 s and an Eddington accretion 
rate. If the infalling material does not possess sufficient an- 
gular momentum, however, a disk will not form, but infall 
of the bound material from the envelope is still likely to 
proceed, albeit in a more spherical fashion. The accretion 
rate during the early phase depends on the details of the 
(yet unclear) explosion mechanism. Estimates by Chevalier 
(1989) yield values in the range 3 x 10" 4 -2 x 10 2 M Q yr" 1 . 
In order for the pulsar mechanism to be able to operate, the 
pressure of the pulsar magnetic field must overcome that 
of the spherical infall. For the accretion rates expected at 
early times, however, the pressure of the accreting material 
dominates over the pulsar pressure even at the neutron star 
surface. Chevalier (1989) estimates that, for accretion rates 
M > 3 x 1O _4 M0 yr -1 , the photon luminosity is trapped by 
the inflow and the effects of a central neutron star are hid- 
den. Once the accretion rate drops below that value, pho- 
tons begin to diffuse out from the shocked envelope; from 
that point on, the accretion rate drops rapidly, and the pul- 
sar mechanism can turn on. Chevalier (1989) estimates that 
this occurs at an age of about 7 months. Therefore, even if 
fallback plays a major role in the initial phase of the SN and 
NS lives, its effects are not expected to be relevant at the 
timescales of interest for the conclusions of this work. 

For the purpose of our analysis, we are especially inter- 
ested in the X-ray luminosity at times long enough so that 
accretion is unimportant, but short enough that rotation is 
still the main source of energy. During a Crab-like phase, rel- 
ativistic particles accelerated in the pulsar magnetosphere 
are fed to a synchrotron emitting nebula, the emission of 
which is characterized by a powerlaw spectrum. Another im- 
portant contribution is the pulsed X-ray luminosity (about 
10% of the total in the case of the Crab) originating directly 
from the pulsar magnetosphere. It should be noted that one 
important assumption of our analysis is that all (or at least 
the greatest majority) of neutron stars goes through an early 
time phase during which their magnetosphere is active and 
converts a fraction of the rotational energy into X-rays. How- 
ever, there is observational evidence that there are objects, 
known as Central Compact Object:^ (CCOs), for which no 
pulsar wind nebulae are detected. Since no pulsations are 
detected for these stars, it is possible that they are simply 
objects born slowly rotating and which hence have a low 
value of Erot- In this case, they would not affect any of our 
considerations, since the L x — E TO t correlation appears to 
hold all the way down to the lowest measured values of L x 
and E ro t. However, if the CCOs are NSs with a high E TO t, 
but for which there exists some new physical mechanism 
that suppresses the magnetospheric activity (and hence the 
X-ray luminosity) to values much below what allowed by 
the scatter in the pulsar L x — E Iot relation, then these stars 
would affect the limits that we derive. Since at this stage 
their nature is uncertain, we treat the all sample of NSs 
on the same footing, although keeping this in mind as a 



9 Examples are central source in Cas A and in Puppis A (e.g. 
Petre et al. 1996; Pavlov et al. 2000). 
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caveat should future work demonstrate the different intrin- 
sic nature of the CCOs with respect to the conversion of 
rotational energy into X-ray luminosity. 

As discussed in §1, for all the neutron stars for which 
both the rotational energy loss, E Iot , and the X-ray luminos- 
ity, L x , have been measured, there appears to be a correla- 
tion between these two quantities. This correlation appears 
to hold over a wide range of rotational energy losses, in- 
cluding different emission mechanisms of the pulsar. Since 
in the high L x regime (young pulsars) of interest here the 
X-ray luminosity is dominated by rotational energy losses, 
the most appropriate energy band for our study is above ~ 2 
keV, where the contribution of surface emission due to the 
internal energy of the star is small. The correlation between 
L x and E to t in the 2-10 keV band was first examined by 
Saito et al. (1997) for a small sample of pulsars, and a more 
comprehensive investigation with the largest sample up to 
date was later performed by P02. They found, for a sample 
of 39 pulsars, that the best fit is described by the relation 

logl^p-io] = 1-34 log£ rot - 15.34 , (1) 

with la uncertainty intervals on the parameters a = 1.34 
and b = 15.34 given by a a — 0.03 and <jj, = 1.11, respec- 
tively. A similar analysis on a subsample of 30 pulsars with 
ages t < 10 6 yr by Guseinov et al. (2004) yielded a best 
fit with parameters a = 1.56, b — 23 A, and correspond- 
ing uncertainties a a = 0.12 and <7b = 4.44. The slope of 
this latter fit is a bit steeper than that of P02; as a re- 
sult, the model by Guseinov et al. predicts a larger fraction 
of high luminosity pulsars from the population of fast ro- 
tating young stars with respect to the best fit of P02. In 
order to be on the conservative side for the predicted num- 
ber of high-L^ pulsars, we will use as our working model 
the one by P02. It is interesting to note, however, that both 
groups find that the efficiency r\ x = L x /E TOt is an increasing 
function of the rotational energy loss _E rot of the star. Fur- 
thermore, the analysis by Guseinov et al. shows that, for a 
given E Iot , pulsars with larger B field have a systematically 
larger efficiency r\ x of conversion of rotational energy into 
X-rays. An increase of r/ x with E TOt was found also in the 
investigation by Cheng, Taam & Wang (2004). They con- 
sidered a sample of 23 pulsars and studied the trend with 
E TO t of the pulsed and non-pulsed components separately. 
Their best-fit yielded Lg ul oc E^ ±0 - 08 for the pulsed com- 
ponent, and L x pul oc E^^ ' 1 for the non-pulsed one. They 
noticed how the former is consistent with the theoretical 
X-ray magnetospheric emission model by Cheng & Zhang 
(1999), while the latter is consistent with a PWN model in 
which L" pul oc -Ef/ t 2 , where p ~ 2 — 3 is the powerlaw in- 
dex of the electron energy distribution. Their best fit for the 
total X-ray luminosity (pulsed plus unpulsed components) 
yielded L x oc £^t , fully consistent with the best-fit 
slope of P02. Along similar lines, recently Li et al. (2007) 
presented another statistical study in which, using Chandra 
and XMM data of galactic sources, they were able to resolve 
the component of the X-ray luminosity due to the pulsar 
from that due to the PWN. Their results were very similar 
to those of Cheng et al. (2004), with a best fit for the pul- 
sar component L X BI oc E^ 0A , and a best fit for the PWN 
(representing the unpulsed contribution) L^ WN oc El o ^ ±0 2 . 
They found that the main contribution to the total luminos- 
ity generally comes from the unpulsed PWN, hence yielding 



the steepening of the Lx — E ro t correlation with E Iot , con- 
sistently with the P02 relation, where the contribution from 
the pulsar and the PWN are not distinguished. For our pur- 
poses, we consider the sum of both contributions, since we 
cannot resolve the two components in the observed sample 
of historical SNe0 

It should be noted that, despite the general agreement 
among the various studies on the trend of L x with E rot , and 
the support from theory that the correlation is expected to 
steepen with E Iot , there must be a point of saturation in or- 
der to always satisfy the condition L x < E Iot . While in our 
simulations we impose the extra condition that rfe < 1, it is 
clear that, until the correlation can be calibrated through 
direct measurements of objects with high values of -Erot, 
there remains an uncertainty on how precisely the satura- 
tion occurs, and this uncertainty is unavoidably reflected in 
the precise details of our predictions. However, unless there 
is, for some reason, a point of turnover above the observed 
range where the efficiency of conversion of rotational energy 
into X-rays turns back into r\ x <C 1, then our general conclu- 
sions can be considered robust. In our analysis, in order to 
quantify the uncertainty associated with the above, we will 
also explore the consequences of a break in r\ x just above the 
observed range (to be mostly conservative). 

Another point to note is that one implicit assumption 
that we make in applying the L x — E Iot relation to very 
young objects is that the synchrotron cooling time t sy nch 
in X-rays is much smaller than the age of the source, so 
that the X-ray luminosity is essentially an instantaneous 
tracer of E TO t- In order to check the validity of this as- 
sumption, we have made some rough estimates based on 
measurements in known sources. For example, let's consider 
the case of the PWN in SN 1986J. Although the field in 
the PWN has not been directly measured, we can use ra- 
dio equipartition and scale it from that of the Crab Neb- 
ula. The Crab's radio synchrotron emission has a minimum 
energy of ~ 6 x 10 48 ergs (see e.g. Chevalier 2005), and 
a volume of ~ 5 x 10 56 cm 3 . The average magnetic field 
is then ~ 550 fiG. We can then scale to the PWN in SN 
1986J, using the fact that the radio luminosity is related to 
that of the Crab by Z/ r ,i986J ~ 200I/ rj crab, while its size is 
about 0.01 times that of the Crab. According to equipar- 
tition, B mln ~ (size) _6/7 L 2/7 (e.g. Willott et al. 1999), so 
this very crude approach suggests that the magnetic field 
in the PWN of SN 1986J is B min ~ 235B C rab ~ 120 mG. 
This yields a very short cooling time in X-rays, i sy iich ~ 5 
hr (assuming a Lorentz factor for the electrons of ~ 10 ), so 
that, if we scale from the Crab nebula, the use of L x — E Iot 
at early times appears reasonable. If, on the other hand, ini- 
tial periods are generally slower than for the Crab, then the 
equipartition energy could be much smaller and the corre- 
sponding lifetimes much longer. Let's then consider a 10 12 G 
pulsar with an initial period of 60 ms (the pulsar produced in 
SN 386 is such a source ). We then have E vot ~ 3 x 10 36 erg/s, 
so that (ignoring expansion losses), the energy deposited in 
the PWN over 20 years would be ~ 2 x 10 45 ergs. For a 
volume similar to that for SN 1986J above, the equiparti- 



10 For typical distances > a few Mpc, the ACIS spatial resolution 
is > 20 pc, while PWN sizes are on the order of a fraction of a 
pc to a few pc. 
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tion magnetic field would be ~ 10 mG, corresponding to 
a lifetime at 2 keV of about 10 days. This is still a short 
enough lifetime for our purposes. Alternatively, for Po ~ 5 
ms and B ~ 10 12 G, we have E vot ~ 6 x 10 40 erg/s, and 
over 20 years, this yields E to t ~ 4 x 10 49 ergs. This im- 
plies B > 1 G, so that tsynch ~ 15 minutes. Therefore, we 
conclude that, overall, the magnetic fields in young PWNe 
are likely strong enough to justify the use of the L x — E Iot 
relation even for the youngest objects in our sample. 

One further point to notice with respect to the L x — E Iot 
correlation is the fact that it is based on a diversity of ob- 
jects. The low end range of the relation, in particular, is 
populated with Millisecond Pulsars (MSPs) , which are spun 
up neutron stars. It is possible that this class of objects 
might bias the correlation of the youngest, isolated pulsars 
in the sample. Generally speaking, once they are spun up, 
the MSPs form PWNe again (e.g. Cheng et al. 2006), and 
the conversion of E Iot into L x , which is practically an in- 
stantaneous relationship (as compared to the ages under 
consideration), should not be dependent on the history of 
the system. The magnetic field of the objects (lower for the 
MSPs than for the young, isolated pulsars), however, might 
influence the conversion efficiency (Guseinov et al. 2004), 
hence biasing the overall slope of the correlation. Overall, in 
our analysis a steeper slope would lead to tighter limits on 
the NS spin birth distribution, and viceversa for a shallower 
slope. What would be affected the most by a slope change 
is the high E TOt tail of the population. Hence, in §4, besides 
deriving results using the L x — E rot for all the pulsars, wc 
will also examine the effects of a change of slope for the 
fastest pulsars. 

The rotational energy loss of the star, under the as- 
sumption that it is dominated by magnetic dipole losses, is 
given by 

b 2 sm 2 en 4 R 6 



E ro t — 



6c 3 



(2) 



where R is the NS radius, which we take to be 10 km, B 
the NS magnetic field, Q = 2tt/P the star angular velocity, 
and 8 the angle between the magnetic and spin axes. We 
take sm6 — 1 for consistency with what generally assumed 
in pulsar radio studies. With sin# = 1 and a constant B 
field, the spin evolution of the pulsars is simply given by 

1 1/2 



P(t) = 



Po 2 + 



m^R^B 

3/c 3 



(3) 



where / « 10 4j gcm 2 is the moment of inertia of the 
star, and Po is its initial spin period. The X-ray lumi- 
nosity of the pulsar at time t (which traces E Iot ) cor- 
respondingly declines as L x — L Xi o(l + t/to)~ 2 , where 
to = 3Ic 3 Pg/B 2 R 6 (2Tv) 2 ~ 6500yr /45Br 2 2 Pio 6p o,io> hav " 
ing defined hs = 7/(10 45 g cm 2 ), R w = R/(W km), 
Po,io = Po/(10 ms). For t<to the flux does not vary signif- 
icantly. Since the ages isN of the SNe in our sample are all 
< 77 yr, we deduce that, for typical pulsar fields, tsN <C to. 
The luminosities of the pulsars associated with the SNe in 
our sample are therefore expected to be still in the plateau 
region, and thus they directly probe the initial birth param- 
eters, before evolution affects the periods appreciably. 

In order to compute the X-ray luminosity distribution 
of a population of young pulsars, the magnetic fields and 
the initial periods of the pulsars need to be known. As dis- 



cussed in §1, a number of investigations have been made over 
the last few decades in order to infer the birth parameters 
of NSs, and in particular the distribution of initial periods 
and magnetic fields. Here, we begin our study by comparing 
the SN data with the results of a pulsar population calcula- 
tion that assumes one of such distributions, and specifically 
one that makes predictions for birth periods in the millisec- 
ond range. After establishing that the SN data are highly 
inconsistent with such short initial spins, we then general- 
ize our analysis by inverting the problem, and performing a 
parametric study aimed at finding the minimum values of 
the birth periods that result in predicted X-ray luminosities 
consistent with the SN X-ray data. 



4 OBSERVATIONAL CONSTRAINTS ON THE 
PULSAR X-RAY LUMINOSITIES FROM 
COMPARISON WITH HISTORICAL SNE 

As a starting point to constrain pulsar birth parameters, 
we consider the results of one of the most recent and com- 
prehensive radio studies, based on large-scale radio pulsar 
surveys, that is the one carried out by Arzoumanian et al. 
(2002; ACC in the following). They find that, if spin down 
is dominated by dipole radiation losses (i.e. braking index 
equal to 3) , and the magnetic field does not appreciably de- 
cay, the magnetic field strength (taken as Gaussian in log) 
has a mean (log Bo [G]) = 12.35 and a standard deviation of 
0.4, while the initial birth period distribution (also taken as 
a log- Gaussian), is found to have a mean (logPo(s)) = —2.3 
with a standard deviation ap > 0.2 (within the searched 
range of 0.1 — 0.7). In the first part of our paper, as a spe- 
cific example of a distribution that predicts a large fraction 
of pulsars to be born with millisecond periods, we use their 
inferred parameters described above. Since in their model 
the standard deviation for the initial period distribution 
is constrained only by the lower limit ap > 0.2, here we 
adopt ap — 0.3. As the width of the velocity dispersion in- 
creases, the predicted X-ray luminosity distribution becomes 
more and more heavily weighed towards higher luminosities 
(sec Figure 1 in Perna & Stella 2004). Therefore, the limits 
that we derive in this section would be even stronger if ap 
were larger than what we assume. We then assume that the 
Lx — E Iot correlation is described by the P02 best fit with 
the corresponding scatter. 

In order to test the resulting theoretical predictions for 
the pulsar distribution of X-ray luminosities against the lim- 
its of the SNe, we perform 10 6 Monte Carlo realizations of 
the compact object remnant population. Each realization is 
made up of iV bj = -/Vsn = 100, with ages equal to the ages 
of the SNe in our sample. The fraction of massive stars that 
leave behind a black hole (BH) has been theoretically esti- 
mated in the simulations by Heger et al. (2003). For a solar 
metallicity and a Salpeter IMF, they find that this fraction 
is about 13% of the total. However we need to remark that 
while, following their predictions, in our Monte Carlo simu- 
lations we assign a 0.13 probability for a remnant to contain 
a BH, the precise BH fraction is, in reality, subject to a cer- 
tain degree of uncertainty. Even taking rigorously the results 
of Heger et al. (2003), one needs to note that their NS vs BH 
fraction (cfr. their fig. 5) was computed assuming a fraction 
of about 17% of Type Ib/c SNe, and 87% of Type II. Our 
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sample, on the other hand, contains about 40% of Type Ib/c 
and 60% of Type II SNe. How the remnant fraction would 
change in this case is difficult to predict. Heger et al. point 
out how normal Type Ib/c SNe are not produced by single 
stars until the metallicity is well above solar. In this case, 
the remnants would be all NSs. At lower metallicities, on 
the other hand, most Type Ib/c SNe are produced in binary 
systems where the binary companion helps in removing the 
hydrogen envelope of the collapsing star. Given these un- 
certainties, while adopting for our simulations the BH/NS 
fraction estimated by Heger et al. for solar metallicity, we 
also discuss how results would vary for different values of 
the BH and NS components. 

If an object is a BH, a low level of X-ray luminos- 
ity (< 10 35 erg/s, i.e. smaller than the lowest measure- 
ment/limit in our SN data set) is assigned to it. This is 
the most conservative assumption that we can make in or- 
der to derive constraints on the luminosity distribution of 
the NS component. If an object is a NS, then its birth 
period and magnetic field is drawn from the ACC distri- 
bution as described above, and it is evolved to its cur- 
rent age (equal to the age of the corresponding SN at the 
time of the observation) with Eq.©. The corresponding X- 
ray luminosity is then drawn from a log-Gaussian distribu- 
tion with mean given by the P02 relation, and dispersion 

<7L X = V Cr a[log£'rot] 2 + of. 

Figure 1 (top panel) shows the predicted distribution of 
the most frequent valuq^j of the pulsar luminosity over all 
the Monte Carlo realizations of the entire sample of Table 1. 
The shaded region indicates the la dispersion in the model. 
This has been determined by computing the most compact 
region containing 68% of the random realizations of the sam- 
ple. Also shown is the distribution of the X-ray luminosity 
(both detections and upper limits) of the SNe (cfr. Table 1). 
Since the measured X-ray luminosity of each object is the 
sum of that of the SN itself and that of the putative pulsar 
embedded in it, for the purpose of this work X-ray detections 
are also treated as upper limits on the pulsar luminosities. 
This is indicated by the arrows in Figure 1. 

Our X-ray analysis, in all those cases where a mea- 
surement was possible, never revealed column densities high 
enough to affect the observed 2-10 keV flux significantly. 
However, if a large fraction of the X-ray luminosity (when 
not due to the pulsar) does not come from the innermost 
region of the remnant, then the inferred Nh would be un- 
derestimated with respect to the total column density to the 
pulsar. The total optical depth to the center of the SN as a 
function of the SN age depends on a number of parameters, 
the most important of which are the ejected mass and its ra- 
dial distribution. The density profile of the gas in the newly 
born SN is determined by the initial stellar structure, as 
modified by the explosion. Numerical simulations of super- 
nova explosions produce density distributions that, during 
the free expansion phase, can be approximated by the func- 
tional form psN = f(v)t~ 3 (see e.g. Chevalier & Fransson 
1994 and references therein). The function f(v) can in turn 



11 For each (binned) value of the pulsar luminosity, we deter- 
mined the corresponding probability distribution resulting from 
the Monte Carlo simulations. The maximum of that distribution 
is what we indicated as the "most frequent" value for each bin. 
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Figure 1. The dashed line shows the distribution of 2-10 keV 
luminosities (either measurements or upper limits) for the entire 
sample of 100 SNe analyzed (upper panel), for the subsample of 
SNe with ages > 10 yr (middle panel), and with ages > 30 yr 
(lower panel). The measured SN luminosities are also treated as 
upper limits on the luminosities of the embedded pulsars. The 
solid line shows the prediction for the X-ray luminosity distri- 
bution of pulsars born inside those SNe, according to the ACC 
birth parameters and the Lx — E TO t P02 relation. The shaded re- 
gions indicate the l-cr confidence level of the model, derived from 
10 6 random realizations of the sample. Independently of the SN 
sample considered, the pulsar luminosity distribution is highly 
inconsistent with the corres^^^^^^t^ 
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be represented by a power-law in velocity, f(v) oc v~™. To 
date, the best studied case is that of SN 1987A. Modeling by 
Arnett (1988) and Shigeyama & Nomoto (1990) yield an al- 
most flat inner powerlaw region, surrounded by a very steep 
outer powerlaw profile, n ~ 9 — 10. For normal abundances 
and at energies belowH 10 keV, Chevalier & Fransson (1994) 
estimate the optical depth at energy Eio = £7/(10 keV) to 
the center of a supernova with a flat inner density profile 
to be r = r s _E ] " 8/ ' 3 Eg^'s! Mf/io*yr 2 > wnere -Esn.si is the 
supernova energy in units of 10 51 erg, and M e j,io is the 
mass of the ejecta in units of 10 M©. The constant t s is 
found to be 5.2 for a density profile with n = 7 in the outer 
parts, and 4.7 for n = 12. From these simple estimates, it 
can be seen that the SN would have to wait a decade or 
so before starting to become optically thin at the energies 
of interest. These estimates however do not account for the 
fact that, if the SN harbors an energetic pulsar in its cen- 
ter, the pulsar itself will ionize a substantial fraction of the 
surrounding neutral material. Calculations of the ionization 
front of a pulsar in the interior of a young SN were per- 
formed by Chevalier & Fransson (1992). In the case of a 
flat density profile in the inner region, and an outer density 
profile with powerlaw n — 9, they estimate that the ioniza- 
tion front reaches the edge of the constant density region 
after a time t yl = 10 tofr E7Jt,4i M lj% s sn,51» where 
_E rot ,4i = iJrot/10 41 erg s _1 , and fi is the fraction of the to- 
tal rotational power that is converted in the form of ionizing 
radiation with a mean free path that is small compared to 
the supernova size. The constant to depends of the composi- 
tion of the core. For a hydrogen-dominated core, to = 1.64, 
for a helium-dominated core, to = 0.69, and for an oxygen- 
dominated core to = 0.28. Once the ionization front has 
reached the edge of the constant density region, the steep 
outer power-law part of the density profile is rapidly ionized. 
Therefore, depending on the composition and total mass of 
the ejecta, an energetic pulsar can ionize the entire mass of 
the ejecta on a timescale between a few years and a few tens 
of years. This would clearly reduce the optical depth to the 
center of the remnant estimated above. 

Given these considerations, in order to make predictions 
that are not as likely to be affected by opacity effects, we 
also performed a Monte Carlo simulation of the compact 
remnant population for all the SNe with ages t > 10 yr, and 
another for all the SNe with ages t > 30 yr. Since the opacity 
scales as t~ 2 , these subsets of objects are expected to be sub- 
stantially less affected by high optical depths to their inner 
regions. The subsample of SNe with ages t > 10 yr contains 
40 objects, while the subsample with ages t > 30 yr contains 
21 SNe. The corresponding luminosity distributions (both 
measurements and limits) are shown in Figure 1 (middle and 
bottom panel respectively), together with the predictions of 
the adopted model (ACC initial period distribution and P02 
L x — E TOt correlation) for the luminosities of the pulsars as- 
sociated with those SN samples. Given the uncertainties in 
the early-time optical depth, we consider the constraints de- 
rived from these subsamples (and especially the one with 
t > 30 yr) more reliable. Furthermore, even independently 
of optical depth effects that can bias the youngest mem- 

12 Above 10 keV, the opacity is dominated by electron scattering, 
which is energy independent. 



bers of the total sample, the subsamples of older SNe have 
on average lower luminosities, hence making the constraints 
on the model predictions more stringent. In the following, 
when generalizing our study to derive limits on the allowed 
initial period distribution, we will use for our analysis only 
the subsets of older SNe. 

In all three panels of Figure 1, the low luminosity tail of 
the simulation, accounting for ~ 15% of the population, is 
dominated by the fraction of SNe whose compact remnants 
are black holes, and for which we have assumed a luminosity 
lower than the lowest SN measurement/limit (~ 10 35 erg/s). 
While it is possible that newly born BHs could be accreting 
from a fallback disk and hence have luminosities as high as a 
few x 10 38 erg/s, our assumption of low luminosity for them 
is the most conservative one for the analysis that we are 
performing, in that it allows us to derive the most stringent 
limits on the luminosity of the remaining remnant popula- 
tion of neutron stars. For these, the high-luminosity tail is 
dominated by the fastest pulsars, those born with periods 
of a few ms. The magnetic fields, on the other hand, are in 
the bulk range of 10 12 - 10 13 G. The law-B field tail pro- 
duces lower luminosities at birth, while the high-B field tail 
will cause the pulsars to slow down on a timescale smaller 
than the typical ages of the SNe in the sample. Therefore, it 
is essentially the initial periods which play a crucial role in 
determining the extent of the high-luminosity tail of the dis- 
tribution. With the birth parameter distribution used here, 
we find that, out of the 10 6 Monte Carlo realizations of the 
sample (for each of the three cases of Fig.l), none of them 
predicts pulsar luminosities compatible with the SN X-ray 
limitsE3 

These results point in the direction of initial periods 
of the pulsar population to be slower than the ms periods 
derived from some population synthesis studies in the radio. 
A number of other investigations in the last few years, based 
on different methods of analysis of the radio sample with 
respect to ACC, have indeed come up to conclusions similar 
to ours. The population synthesis studies of Faucher-Giguere 
& Kaspi (2005) yielded a good fit to the data with the birth 
period described by a Gaussian with a mean period of 0.3 s 
and a spread of 0.15 sec. Similarly, the analysis by Ferrario 
& Wickramasinghe (2006) yielded a mean period of 0.23 sec 
for a magnetic field of 10 12 G. We performed Monte Carlo 
simulations of the X-ray pulsar population using the birth 
parameters derived in those studies above, and found them 
to be consistent with the SNe X-ray limits shown in Figure 1. 

In order to generalize our analysis beyond the testing of 
known distributions, we performed a number of Monte Carlo 
simulations with different initial spin period distributions 
and a mean magnetic field given by the optimal model of 
Faucher-Giguere & Kaspi (2006). This is a log-Gaussian with 
mean (log(S/G)) = 12.65 and dispersion <Ti og s = 0.550. 

13 We need to point out that, in the study presented here, we 
refrain from performing detailed probability analysis. This is be- 
cause, given the observational uncertainties of some of the input 
elements needed for our study (as discussed both above and in 
the following), precise probability numbers would not be espe- 
cially meaningful at this stage. 

14 The inferred values of the magnetic field in different studies 
are all generally in this range, even for very different inferred 
spin birth parameters. Furthermore, Ferrario & Wickramasinghe 
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Figure 2. Fraction f c of Monte Carlo realizations of the SN 
sample for which the 2-10 keV luminosities of the pulsars are be- 
low the limits of the corresponding SNe. This is shown for different 
distributions of the initial spin periods, described by Gaussians 
of mean Po and dispersion <jp Q . In the upper panel, the sample 
includes only the SNe of ages > 10 yr (cfr. Fig.l, middle panel), 
while in the lower panel only the SNe with ages > 30 yr (cfr. 
Fig.l, lower panel) are included for the montecarlo simulations. 
Independently of the sample considered, initial periods Po % 40 
ms are inconsistent with the SN data. 

Figure 2 shows the fraction fg im of the montecarlo sim- 
ulations of the SN sample for which the luminosity of each 



(2006) note that the pulsar birth period that they infer is almost 
independent of the field value in the range log B(G) = 10 — 13, 
where the vast majority of isolated radio pulsar lie. 



pulsar is found below that of the corresponding SN. The 
sample of SNe selected is either the one with ages > 10 
yr (top panel), or the one with ages > 30 yr (bottom 
panel), which allow tighter constraints while minimizing op- 
tical depth effects. Monte Carlo realizations of the samples 
have been run for 50 Gaussian distributions of the period 
with mean in the range 20 — 100 ms, and, for each of them, 
4 values of the dispersion^! op Q between 1 and 50 ms. For 
each value of the period, we performed 100,000 random re- 
alizations0. Details of the results vary a bit between the 
two age-limited subsamples. This is not surprising since the 
extent to which we can draw limits on the pulsar periods de- 
pends on the measurements/limits of the X-ray luminosities 
of the SNe in the sample. In a large fraction of the cases, we 
have only upper limits, and therefore our analysis is depen- 
dent on the sensitivity at which each object has been ob- 
served. Independently of the sample, however, we find that, 
for initial periods Prj ^ 35 — 40 ms, the distribution of pulsar 
luminosities is highly inconsistent with the SN data for any 
value of the assumed width of the period distribution. 

We need to emphasize that the specific value of /° im 
as a function of Po should be taken as representative. Vari- 
ous authors have come up with slightly different fits for the 
L x — E Iot correlation. If, for example, instead of the fit by 
Possenti et al. (2002) we had used the fit derived by Gu- 
seinov et al. (2004), then the limits on the period would 
have been more stringent. On the other hand, if, for some 
reason, the efficiency rj x of conversion of E TOt into L x be- 
comes low at high E Iot , then our results would be less con- 
straining. In order to assess the robustness of our results 
with respect to changes in rj x , we ran simulations of the pul- 
sar population assuming that, for P rot > E™*- ohs ~ 10 39 
erg s _1 (where P™j* x ' obs is the maximum observed E rot ), the 
efficiency becomes r( x = erj x , and we tried with a range of 
values of e < 1. This test also addresses the issue of a bias 
in our results deriving from a possible shallower slope for 
the youngest pulsars of the population, as discussed in §3. 
We ran Monte Carlo simulations for the ACC birth parame- 
ters, and decreased e by increments of 0.02. We found that, 
only for the very low value e ~ 10 -4 , a sizable fraction of 
~ 5% of the simulations predicts pulsar X-ray luminosities 
that are fully consistent with the SN data. Therefore, we 
conclude that our results on the millisecond birth periods of 
pulsars are reasonably robust with respect to uncertainty in 
the L x — E TO t for the youngest members of the population. 

Another systematic that might in principle affect our 
results would arise if a fraction of neutron stars is born with 
a non-active magnetosphere so that their X-ray luminosity 
at high P rot is much smaller than for the active pulsars, 



5 The dependences with ap Q should be taken as representative 
of the general trend, since it is likely that ap Q and Po might be 
correlated. But since no such correlations have been studied and 
reported, we took as illustrative example the simplest case of a 
constant op for a range of Po- 

16 The number of random realizations is smaller here with respect 
to Fig.l for computational reasons since, while each panel of Fig.l 
displays a Monte Carlo simulation for one set of parameters only, 
each panel of Fig. 2 is the results of 200 different Monte Carlo 
realizations. For a few cases, however, we verified that the results 
were statistically consistent with those obtained with a larger 
number of random realizations. 



© 0000 RAS, MNRAS 000, 000-000 



How rapidly do neutron stars spin at birth? 11 



then the limits on the initial periods of the "active" pulsars 
would be less stringent. An example of non-active neutron 
stars could be that of the CCOs discussed in §3. However, 
until the fraction of these stars becomes well constrained by 
the observations and an independent L x — E Iot is established 
for them, it is not possible to include them quantitatively 
in our population studies. Similarly, the precise fraction of 
BHs versus NSs in the remnant population plays a role in 
our results. A larger fraction of BHs would alleviate our con- 
straints on the initial spin periods, while a smaller fraction 
would, obviously, make them tighter. If a fraction of those 
BHs had a luminosity larger than the maximum assumed up- 
per limit in our simulations (due to e.g. accretion from a fall- 
back disk as discussed above), then our results would again 
be more constraining. While our work is the first of its kind 
in performing the type of analysis that we present, future 
studies will be able to improve upon our results, once the 
possible systematics discussed above are better constrained, 
and deeper limits are available for the full SN sample. 



5 SUMMARY 

In this paper we have proposed a new method for prob- 
ing the birth parameters and energetics of young neutron 
stars. The idea is simply based on the fact that the X-ray 
measurements of young supernovae provide upper limits to 
the luminosity of the young pulsars embedded in them. The 
pulsar X-ray luminosity on the other hand, being directly re- 
lated to its energy loss, provides a direct probe of the pulsar 
spin and magnetic field. Whereas pulsar birth parameters 
are generally inferred through the properties of the radio 
population, the X-ray properties of the youngest members 
of the population provide tight and independent constraints 
on those birth parameters, and, as we discussed, especially 
on the spins. 

The statistical comparison between theoretical predic- 
tions and the distribution of X-ray luminosity limits that we 
have performed in this work has demonstrated that the two 
are highly inconsistent if the bulk of pulsars is born with pe- 
riods in the millisecond range. Whereas we cannot exclude 
that the efficiency r\ x of conversion of rotational energy into 
X-ray luminosity could have a turnover and drop at high 
values of E Iot to become rj x <C 1, the 2-10 keV pulsar data 
in the currently observed range of E rot do not point in this 
direction (but rather point to an increase of r) x with E TO t), 
and there is no theoretical reason for hypothesizing such a 
turn over. However, even if such a turnover were to exist 
just above the observed range of E TOt , we found that only 
by taking an efficiency r)' x ~ 10~ 4 r] x above E^: x '° bs , our 
results would lose their constraining value for the ms spin 
birth distributions. Therefore, we can robustly interpret our 
results as an indication that there must be a sizable fraction 
of pulsars born with spin periods slower than what has been 
derived by a number of radio population studies as well as 
by hydrodynamic simulations of SN core-collapse (e.g. Ott 
et al 2006). Our findings go along the lines of a few di- 
rect measurements of initial periods of pulsars in SNRs (see 
e.g. Table 2 in Migliazzo et al. 2002), as well as some other 
population synthesis studies (Faucher-Gigure & Kaspi 2006; 
Ferrario & Wickramasinghe 2006; Lorimer et al. 2006). Our 
results for the bulk of the pulsar population, however, do not 



exclude that the subpopulation of magnetars could be born 
with very fast spins, as needed in order to create the dynamo 
action responsible for the B-field amplification required in 
these objects (Thompson & Duncan 1993). Because of their 
very short spin-down times, the energy output of magnetars 
can be dominated by the spin down luminosity only up to 
timescalcs of a fraction of year, during which the SN is still 
too optically thick to let the pulsar luminosity go through. 
Therefore, our analysis cannot place meaningful constraints 
on this class of objects. 

Finally, our results also bear implications on the con- 
tribution of young pulsars to the population of the Ultra 
Luminous X-ray sources (ULXs) observed in nearby galax- 
ies. The model in §3 predicts that a sizable fraction of that 
population could indeed be made up of young, fast rotat- 
ing pulsars (Perna & Stella 2004). However, the analysis 
performed here shows that the contribution from this com- 
ponent, although expected from the tail of the population, 
cannot be as large as current models predict. 

The extent to which we could perform our current anal- 
ysis has been limited by the size of the SNe sample, and, es- 
pecially, by the available X-ray measurements. The fact that, 
in a large fraction of the sample, we have limits rather than 
detections, means that a large improvement can be made 
with deeper limits from longer observations. The deeper the 
limits, the tighter the constraints that can be derived on the 
spin period distribution of the pulsars. The analysis pro- 
posed and performed here is completely uncorrelated from 
what done in radio studies, and therefore it provides an 
independent and complementary probe of the pulsar spin 
distribution at birth (or shortly thereafter) ; our results pro- 
vide stronger constraints on theoretical models of stellar core 
collapse and early neutron star evolution, making it even 
more necessary to explain why neutron stars spin down so 
rapidly immediately after birth (see also Thompson, Chang 
& Quataert 2004; Metzger, Thompson & Quataert 2007). 
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